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Abstract 

We compute the contributions to continuum photon pair production at hadron colhders from 
processes initiated by gluon-gluon and gluon-quark scattering into two photons through a four-leg 
virtual quark loop. Complete two-loop cross sections in perturbative quantum chromodynamics 
are combined with contributions from soft parton radiation resummed to all orders in the strong 
coupling strength. The structure of the resummed cross section is examined in detail, including 
a new type of unintegrated parton distribution function affecting azimuthal angle distributions 
of photons in the pair's rest frame. As a result of this analysis, we predict diphoton transverse 
momentum distributions in gluon-gluon scattering in wide ranges of kinematic parameters at the 
Fermilab Tevatron and the CERN Large Hadron Collider. 
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I. INTRODUCTION 



Advances in the computation of higher-order radiative contributions in perturbative quan- 
tum chromodynamics (PQCD) open opportunities to predict hadronic observables at an 
unprecedented level of precision. Full realization of this potential requires concurrent im- 
provements in the methods for QCD factorization and resummation of logarithmic enhance- 
ments in hadronic cross sections in infrared kinematic regions. All-orders resummation of 
logarithmic corrections, such as the resummation of transverse momentum (Qt) logarithms 
in Drell-Yan-like processes [l], is increasingly challenging in multi-loop calculations as a 
result of algebraic complexity and new types of logarithmic singularities associated with 
multi-particle matrix elements. 

In this paper, we address new theoretical issues in Qt resummation at two-loop accuracy. 
We focus on photon pair production, particularly on the gluon-gluon subprocess, gg 77, 
one of the important short-distance subprocesses that contribute to the inclusive reactions 
PP ~^ •yyX at the Fermilab Tevatron and pp 77X at the CERN Large Hadron Collider 
(LHC). This hadronic reaction is interesting in its own right, and it is relevant in searches 
for the Higgs boson h, where it constitutes an important QCD background to the pp 
hX —>■ 77X production chain [1, 0, A reliable prediction of the cross section for gg —>■ 77 
is needed for complete estimates of the 77 production cross sections, a task that we pursue 
in accompanying papers [sl, @|. 

The lowest-order contribution to the cross section for gg 77 arises from a 2 ^ 2 
diagram of order involving a 4-vertex virtual quark loop [Fig. [I](a)]. We evaluate 

all next-to-leading (NLO) contributions of order C(a^a^) to the gg 77 process shown 
in Figs. [I](b-e). An important new ingredient in this paper is the inclusion of the gq 
j'-fq process. Fig. [I](d), a necessary component of the resummed NLO contribution. Our 
complete treatment of the NLO cross section represents an improvement over our original 
publication in which the large-Q^ behavior of the gg subprocess was approximated, and 
the gq contribution was not included. Furthermore, we resum to next-to-next-to-leading 
logarithmic (NNLL) accuracy the large logarithmic terms of the form ln(Q|,/Q^) in the 
limit when Qt of the 77 pair is much smaller than its invariant mass Q. Our NNLL cross 
section includes the exact C coefficients of order as for gg + gq ^ llX, and the functions 
A and B of orders and in all subprocesses, with these functions defined in Sec. [Tll 

We begin in Sec. II with a summary of kinematics and our notation, and we outline the 
partonic subprocesses that contribute to 77 production. In this section, we also derive a 
matrix element for the qg 77(7 process shown in Fig.[l](d), a subprocess whose contribution 
is required to obtain consistent resummed predictions for all values of Qt- We obtain the 
C(a^a^) cross section for the gq •-y^q process from the color-decomposed qqggg amplitudes 
in Ref. |8|. 

The rich helicity structure of the gg 77 matrix element is addressed in Sec. III. The 
helicity dependence requires a new type of transverse-momentum dependent (TMD) parton 
distribution function (PDF) associated with the interference of amplitudes for initial-state 
gluons of opposite helicities. The existence of the helicity-flip TMD PDF modifies the 
azimuthal angle distributions of the final-state photons, an effect that could potentially 
be observed experimentally. By contrast, in vector boson production pp ' VX (with 
V = 7*, VF, Z, ...), such helicity-flip contributions are suppressed as a result of the simple 
spin structure of the lowest-order qqV coupling. In this section, we establish the presence 
of helicity interference in the finite-order 2^3 cross sections by systematically deriving 
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Figure 1: Representative parton scattering subprocesses for diphoton production in gluon-gluon 
scattering. 



their soft and collinear limits in the splitting amplitude formalism 

We show how the helicity-flip TMD PDF arises from the general structure of the small-Qr 
resummed cross section. 

Section IV contains some numerical predictions for the Tevatron and LHC, where we 
show the fraction of the rate for 77 production supplied by the gg + gq subprocess. The 
generally expected prominence of gg + gq scattering at the LHC is only partially supported 
by our findings. The large gg partonic luminosity cannot fully compensate for the small 
cross section associated with gg scattering. Our findings are summarized in Sec. V. Three 
Appendices are included. In Appendix [Al we present some of the details of our derivation 
of the amplitude for the subprocess qg 'J'jg- In Appendix [Bl we derive the small-Qr 
asymptotic form of the NLO cross section for — > 77. 



II. NOTATION AND SUBPROCESSES 



Notation 



We consider the scattering process /ii(Pi) + /i2(-P2) — * 7(-P3) + tI-Ri) + where hi and 
/12 are the initial-state hadrons. In terms of the center-of-mass collision energy v^, the 77 
invariant mass Q, the 77 transverse momentum Qt, and the 77 rapidity y, the momenta 
and P2 of the initial hadrons and q^ = + of the pair are expressed in the laboratory 
frame as 



Pf 



-'2 
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{1,0,0,1}; 
{1,0,0,-1}; 



+ cosh y, Qt, 0,JQ^ + Q\ sinh y 



The light-cone momentum fractions for the boosted 2^2 scattering system are 



(1) 

(2) 
(3) 



2(P2,i ■ q) VQ' + QtC^" 



S 



(4) 



Decay of the 77 pairs is described in the hadronic Collins-Soper frame [15|. The Collins- 
Soper frame is a rest frame of the 77 pair (with q^ = {Q, 0,0,0} in this frame), chosen 
so that (a) the momenta Pi and P2 of the initial hadrons lie in the Oxz plane (with zero 
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azimuthal angle), and (b) the z axis bisects the angle between Pi and — ^2- The photon 
momenta are antiparallel in the Collins-Soper frame: 

= y {0, sin0* cosv9^,, sin^* sinv?^,, cos^^*} , (5) 

Pi = y {0, — sin6'* cos(y9*, — sin^* sin(y9^,, — cos6'*} , (6) 

where 6^, and (p^, are the photon's polar and azimuthal angles. Our aim is to derive re- 
summed predictions for the fully differential 77 cross section da / {dQ'^dydQ'^dVL^), where 
dVL^, = d cos O^dip^ is a solid angle element around the direction of P3 in the Collins-Soper 
frame of reference defined in Eq. (I5|). The parton momenta and helicities are denoted by 
lowercase pi and Aj. 



B. Scattering contributions 

We concentrate on direct production of isolated photons in hard QCD scattering, the 
dominant production process at hadron colliders. A number of hard-scattering contributions 
to the processes qq + g o — » 77, as well as photon production via fragmentation, have been 



studied in the past Il8l |. Our numerical calculations include the lowest-order process 



qq 77 of order C(a^) and contributions from qq — > 77(7 and g's' ~^ TtV of order O(a^as), 
where = e^/in and ^^(/i) = g^/4:'7i are the running QED and QCD coupling strengths. 

Glue-glue scattering is the next leading direct production channel, with the full set of 
NLO contributions shown in Fig. [H Production of 77 pairs via a box diagram in gg scat- 
tering as in Fig. [I](a) is suppressed by two powers of as compared to the lowest-order 
qq — > 77 contribution, but is enhanced by a product of two large gluon PDF's if typical 
momentum fractions x are small. The main O(a^a^), or NLO, corrections, include one-loop 



gg ^^g diagrams (b) and (c) derived in [20, mf, as well as 4-leg two-loop diagrams (e) 
computed in 23|. The real and virtual diagrams are combined in Ref. to obtain the full 



NLO contribution from gg scattering. In this study we also include subleading NLO contri- 
butions from the process (d), gqs 'J'jqs via the quark loop, where qs = J2i=u,d,s,...i'li + Qi) 
denotes the fiavor-singlet combination of quark scattering channels. The gqs — ^ Jiqs helicity 
amplitude is derived from the one-loop qqggg amplitude and explicitly presented in Ap- 
pendix [Al As a cross check, we verified that this amplitude correctly reproduces the known 
collinear limits. Our result does not confirm an expression for this amplitude available in 
the literature [Hi, which does not satisfy these limits. When evaluated in our resummation 
calculation under typical event selection conditions, gg + gqs scattering contributes about 
20% and 10% of the total rate at the LHC and the Tevatron, respectively, but this fraction 
can be larger in specific regions of phase space. 



III. THEORETICAL PRESENTATION 



A. Small-Qr asymptotics of the next-to-leading order cross section 

When the transverse momentum Qt of the diphoton approaches zero, the NLO produc- 
tion cross section da/ {dQ^dy dQ'^dfl^,), or briefiy P{Q,QT,y,^*), is dominated by 77 recoil 
against soft and collinear QCD radiation. In this subsection we concentrate on the effects 
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of initial-state QCD radiation and derive the leading small-Qr part of the NLO differential 
cross section, called the asymptotic term A{Q,QT,y, fi*). 

The 0{as) asymptotic cross section valid at Q"^ <^ consists of a few generalized 
functions that are integrable on an interval < Qt < Pt, with Pt being a finite value of 
transverse momentum: 



AiQ, Qt, y, fi*) = F,(Q, y, fi,)(5(QT) 



+ F,{Q,y,n, 



4ln£ 



T 



+ Fo{Q,y,n,) 



Q 



T. 



+ .... 



(7) 



The "+" prescription [fiQr)], is defined for a function fiQr) and a smooth function g^Qr) 



as 



p2 



dQl [f{QT)].9{QT) 



dQlfiQr) {9{QT)-9m 



[/(Qt)] , = /(Qt) for Qt + 0. 



(8) 
(9) 



Subleading terms proportional to {Q /QtY with p < 1 are neglected in Eq. ((71). Its form is 
infiuenced by spin correlations between the initial-state partons and final-state photons. As 
a consequence of these spin correlations, the functions F^, Fq, and Fi depend on the direction 
of the final-state photons in the Collins-Soper frame (the polar angle 6** and sometimes the 
azimuthal angle y?*). 

The spin dependence of the small-Qr cross section in the gg 77(7 and gqs 
^^qs channels is complex. The Born-level process g{pi,Xi) + (7(^2,^2) 7(^3,^3) + 
7(p4, A4) is described by 16 non-zero helicity amplitudes Al4(pi, Ai; ^2, -^2; Ps, A3; p4, A4) = 
Al4(Ai, A2, A3, A4) for quark-box diagrams of the type shown in Fig.[T](a). The normalization 
of J^4{Xi, A2, A3, A4) is chosen so that the unpolarized Born gg 77 cross section reads as 



da, 



gg 



dQ'^dy dQTdVL^ 



S{Qi 



Born 



s 



fg/hi{Xl, flF)fg/h2{x2, fJ'P), 



where 
with 

and 



a, 



(0) 



a\Q)al{Q) 



327rQ2(Ar2 _ 1) 



(10) 

(11) 

(12) 
(13) 



E l-M4(Ai,A2,A3,A4)r. 

Al,A2,A3,A4=±l 

In these equations, A^^c = 3 is the number of QCD colors, is the fractional electric charge 
(in units of the positron charge e) of the quark i circulating in the loop, and fg/h^x^iip) 
is the gluon PDF evaluated at a factorization scale ^ip- The right-hand side of Eq. (fT3l) 
includes summation over gluon and photon helicities Aj, with « = 1, ...,4. 

At NLO, the small-Qr cross section is proportional to the angular function Sc,(^^*) (the 
same as in the Born cross section), and another function 



S;(^„ <^,) = E -^K^i' ^2, A3, A4)A^4(-Ai, A2, A3, A4) 

Ai,A2,A3,A4=±l 

= a(°)L'(^,)cos2yp,. 



(14) 



5 



140 

120 

100 

80 

60 

40 

20 


0.3 

0.2 

0.1 




0.2 



0.4 0.6 

(1+Cos [(9*J)/2 



0.8 



Figure 2: The functions Lg{9^) and L'g{9^) arising in the gg —>■ 77 asymptotic cross section ([T6] 
and their ratio. 



The function S^(6'*,</9*) is obtained by spin-averaging the product of the amplitude 
Aii{Xi, A2, A3, A4), and the complex-conjugate amplitude A4.l{—Xi, A2, A3, A4) evaluated with 
the reverse sign of the helicity Ai. The sign flip for Ai results in dependence of S^(6'^., </9*) 
on cos2(y9*. The 9^ dependence of ^'g{9^, ip^) enters through the function 



-4Re (Mil., + ^ff 1,1,-1 + ^-'1,1,1,-1 + 1) > 



(15) 



presented in terms of reduced amplitudes M^^,^^^^^^^^ in the notation of Ref. 22 1. For com- 
parison, the functions Lg{9^) and L'g{9^) are plotted versus (1 -f- cos6'*)/2 in Fig. [21 

The NLO asymptotic term in the sum of the contributions from the gg 77 and 
gqs — > 77 channels (denoted as gg + gqs channel) is 



A{Q,QT,y,n,) = U^g{9:, 



6{QT)Fg4Q,y,9,) + Fg^+{Q,y,QT] 



(16) 



Here 



Fg,5 = fg/hA^l,I^F)fg/h2iX2,I^F) + 2^h^^\9^ 



+ 



TV 



c 



(l,c) 



g/a W fa/hi\ (Xufip) - [Pg/a ® fa/h^\ Hf) ~Q j /9/?»2 (3^2 , /^f) 



+fg/hA^l, I^f) [\f-^gi^ ® /aMaJ I^f) - [Pg/a ® /a/feaj (^2, I^f) In |; (17) 

1 



lip 



9,H 



^-^{fa/hi {Xl,f^F)fg/h2{x2, f^F) (^3^'^^ 



— In ^ 
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+ 



Q't 



( Pg/a ® fa/hi {Xl, fip) fg/h2{x2, I^f) 
+ fg/hi /If) Pg/a ® fa/h2 {X2, I^f) 



(18) 



and 



F' = — ^ 



" 1 ■ 








( 

+ 





+/a//ii(a;i,/iF) Pg/g®fg/h2 {x2,^^f))- 



(19) 



The C(a,/7r) coefficients A'i^''\ iSj^'") and functions C^);;;^ (x, 6) , h^g \e,) are defined and 

listed explicitly in Ref. The function h^piO^) denotes an 0(Q;s/vr) correction to the hard- 
scattering contribution Ti in the resummed cross section, cf. Sec. IIII D[ The convolutions 



g/a 



fa 



/h 



and 



CgiJ ® fa/h , defined for two functions f{x,jj,p) and g{x,iJ,F) as 
[f ®g]{x,fiF) = 



y/(^,M^(|,/^F), 

are summed over the intermediate parton's fiavors a = g, qs (gluon and the fiavor-singlet 
combination of quark-scattering channels). In addition to the conventional splitting func- 
tions Pg/g{x) and Pg/qg{x) arising in Fg^+, a new splitting function 



where Ca = Nc 



P'^,^{x) = 2Ca{1-x)/x, (20) 
3, enters the (/^^.-dependent part of the asymptotic cross section through 



F' 

For completeness, the small-Qr asymptotic form Eq. (fT6l) for the gg + gqs channels is de- 
rived in Appendix [Bl The existence of the 99 dependent singular contribution proportional 
to S^(6'*, ip^.) is established by examining the factorization of the 2^3 cross section in the 
limit of a collinear gluo n emission. It follows directly from factorization rules for helicity 
amplitudes 0, 0, [13, H, [H, El, as well as from the dipole factorization formalism H . 

In contrast, the NLO quark-antiquark contribution qq — > 77 does not include a spin-fiip 
contribution, as a result of the simple structure of the Born contribution in qq scattering 
(see also Sec. IIII CI) . 



B. Resummation 



To predict the shape of do jdQp distributions, we perform an all-orders summation of 
singularities ^(Qt) and Q^^ln^ (Q^/Qt) ^ in the asymptotic cross section, which coincides 
with the perturbative expansion of the resummed small-Qr cross section obtained within the 
Collins-Soper-Sterman formalism [ll, [26, 27|. In this formalism, we write the fully differential 
cross section as 



da{hih2 77) 
dQ dQ\ dy dfl^. 



WiQ,QT,y, n,) + YiQ,QT,y, fi. 



(21) 
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Figure 3: The structure of the resummed form factor W{Q,b,y,Q^). 



The term W contains large logarithmic contributions of the form In^ {Q/Qt) from initial- 
state radiation, while Y is free of these logs and calculated using collinear QCD factorization 
(cf. the end of Sec. [TTlD| . 

The function W may be expressed as a Fourier-Bessel transform of a function 
W{Q,b,y,Q^) in the impact parameter (6) space, 

^e^«-''W^(Q,6,y,fi,). (22) 

The generic form of W{Q,h,y,Vt^) in the qq + qg — > 77 and gg + gqs — > 77 channels 
can be determined by solving evolution equations for the gauge- and renormalization-group 
invariance of W{Q, b, y, Q^): 

wiQ,b,y,n,) = J2 E H^'^'^'^(Q,^^*)(H^'^''^'^'^(g,a 

a Ai,a;,A2,A^,A3,A4 

x^';^ (xi, b) b) e-^»(«'^). (23) 

It is composed of the hard-scattering function T-C^^^^^'-^^^{Q,Q^:) and its complex conjugate, 

(A' A' A' A4 \ * 

Ti-a^ ^ ^{Q, fi*) j ; the Sudakov exponential exp {—Sa{Q, &)); and parton distribution ma- 
trices Vl]l]{xi,b). 

The multiplicative structure of Eq. ( l23l l reflects the topology of the dominant cut diagrams 
in the small-Qr cross sections shown in Fig. [3l The function 7-^^i^2A3A4 (describes the hard 
2 — i> 2 scattering subprocess a(pi, Ai) +a(p2, A2) liPs, A3) + 7(^4, A4), with a = u,u, d, d, ... 
in qq 77, and a = a = g in gg ^ 77. All momenta in H have virtualities of order Q^. 
For now, we consider the leading contribution to 7-^^i^2A3A4^ which reads as 'j-C^^^iMU _ 



cri°^7V(4(Ai, A2, A3, A4), where the Born helicity amplitude A^4(Ai, A2, A3, A4) and overall 
constant normalization a^^ are introduced in Sec. IIII A[ Sometimes 7Y^i^2A3A4 ^^so includes 
finite parts of higher-order 2 — > 2 virtual corrections, as discussed in Sec. IIII D[ 

Similarly, (ny'^''^\Q,n,)) arises from the complex-conjugate amplitude 
A^KAi, A2, A3, A4) and possible loop corrections to it. The helicities A'^^ and Ag in 
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Ha^ 2 3 4j j^ggj^ j^Q^ coincide with Ai and A2 in 7Y^i-^2A3A4 'jj^g right-hand side of Eq. f l23ll 
is summed over flavors a and helicities A^, A';, of the partons entering HH*, as well as over 
helicities A3 and A4 of the final-state photons. 
The Sudakov exponent 



Sa{Q,b) 



Aa (Ci,/x) In 



(24) 



resums contributions from the initial-state soft and soft-collinear gluon emission (indicated 
by gluon lines connecting e""^ to 7i, 7i*, and Va/hi^ib) in Fig. [3l). Here Ci and C2 are 
constants of order unity. The functions Aa{Ci,Ji) and Ba{Ci,C2, ft) can be evaluated in 
perturbation theory at large scales /i^ ^ ^qcdi hence for large Q and small b. 

The collinear emissions are described by parton distribution matrices V^/'l^{x,b), where 
A and A' denote the helicity state of the intermediate parton a to the left and right of 
the unitarity cut in Fig . [3l The matrix ^^/^(a:, b) is derived from a matrix element of the 



light-cone correlator [28|, l29|, l30|, |31|, |32| for finding parton a inside the parent hadron h. 

It is convenient to introduce sums of diagonal and off-diagonal entries of the helicity 
matrix V^fl^{x,b), 



Va/H{x,b) = Y.V^'ix,b) 



(25) 



and 

K/hi^^b) 

In this notation, Eq. (l23l) can be rewritten as 



(26) 



W{Q,b,y,n, 



'Sa{Q,b) 



S 



j:a{6*)Va/hA^l,b)Va/h2{x2,b) + 



+ Ki^*,f*) [^a//M(a;i,&)^a/ft2(a;2,&) +PaMi(a;i,6)Ps/h2(a;2,&) 
+ S"(^,,(/?*)P^/;,^(xi,6)p;/^^(x2,6) \, 



(27) 



where 



and 



-T/Ai A2A3A4 
''■a 



E 

Ai,A2,A3,A4 

'y ^ -XVA1A2A3A4 

Ai,A2,A3,A4 



-T/Ai A2A3A4 / ni 
' '-a \ ' ^1 



-A1A2A3A4 



Ai,A2,A3,A4 



A1A2A3A 



4 (^^-A,- 



A2A3A4 



(28) 
(29) 

(30) 



The unpolarized parton distribution Va/h{x, b) coincides with the Fourier-Bessel transform 
of the unpolarized transverse-momentum-dependent (TMD) parton density Va/h{x, kx) [Hi 
for finding parton a with light-cone momentum fraction x and transverse momentum kx- At 
small 6, Va/h{x, b) is reduced to a convolution of unpolarized fc-r— integrated parton densities 
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fa/h{x,fi) and Wilson coefficient functions Ca/a'{.x,h]C /C2, jj), evaluated at a factorization 
scale /i of order 1/6: 



QCD 



E 



' X Ci 



(31) 



Perturbative entries with Aj = reduce in total to the product of the unpolarized Born 
scattering probability and unpolarized resummed functions: 



WiQ,b,y,n, 



A.=A' 



r, S 



The function ^g{0^) is shown explicitly in Eq. (fTTj) . 



®/c2M2 (a;2,6;/i)- (32) 



C. Spin-flip term in gluon scattering 



We concentrate in this subsection on the spin-flip distribution V'g/f^{x, h) in gluon scatter- 
ing. Its existence is warranted by basic symmetries of helicity- and transverse-momentum- 
dependent gluon distribution functions |3j|. This function, which describes interference of 
the amplitudes for nearly collinear gluons with opposite helicities, coincides with the func- 
tion H-^ in Ref. 3j| up to an overall factor. It contributes to unpolarized Qt distributions, 



(with 



^AiA2A3A4 gj^gj^ ^J^g 

-A'2 . The 



because the hard-scattering product '^^i-^2A3A4 

quark box helicity amplitude in Fig. [T](a)) does not vanish for Ai = — A'^ or A2 
presence of 'Pg/f^{x,b) modifies dependence of the resummed cross section on the photon's 
azimuthal angle v?* in the Collins-Soper frame. It vanishes after the integration over is 
performed. In contrast, the helicity-diagonal part of W{Q, b, y, Q^:) is independent of v?*, cf. 
Eq. 

The gluon function 'Pg^f^{x, b) is invariant under time reversal (i.e., is T-even) and acquires 

large contributions proportional to the unpolarized T-even PDF's Vg/h{x, b) in the process 
of gluon radiation. These contributions require resummation via PDF evolution equations 
(similar to Dokshitzer-Gribov-Lipatov-Altarelli-Parisi equations 35, 36, [13, 11^) in order to 
predict the ip^: dependence in the gg channel. 

At one loop, the mixing of spin-flip and unpolarized gluon PDF's is driven by the convo- 
lution Pg^g ® fg/h {x.hf) of the spin-flip splitting function Pg^g{x) shown in Eq. (I20ll with 
the gluon PDF fg/h{x,np). This convolution may be comparable to or exceed the analogous 
convolution Pg/g ® fg/h (a;, /xf) of the unpolarized splitting function Pg/g{x) for some x and 
Hf values, as shown in Fig. [H As a result of the mixing, an additional (y9=i,-dependent term 



V5*j Ols 



2'kSQ 



P' 



fg/hi {Xl, fip) fg/h2{X2, fJ'P) 

+fg/hi (a;i,/iF 



P'g/g ® fg/h^ 



X2,^lF, 



(33) 



arises in the unpolarized 0{as) asymptotic piece, cf. Eq. (ITGll . It is produced by 
the perturbative expansion of the entry proportional to Tj'g{9tf^(pj,)Vg/h{xi,b)V'gij^{xj^b) in 
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W{Q,b,y,Q^), with S^(0*,(y9*) shown explicitly in Eqs. (fT4ll and (fTSll . Generally, the 
dependent contribution is not small, even though it is suppressed comparatively to the 
unpolarized collinear contribution by the ratio L'g{6^) / Lg{6^) shown in Fig. [2l For example, 
for Q = 100 GeV at the LHC, its magnitude constitutes up to about a half of the collinear 
unpolarized asymptotic contribution, 



a. 



([^g/a ® /a/hi] {Xi, jJ^p) fg/h2{x2, /J'f) + fg/hi{Xl, ^f) [Pg/a ® fa/hi] {x2,fJ'F))- (34) 



27lSQ 

The 0{as) spin-flip gg contribution does not mix with the gqs contribution. 

deflned in H 



In terms of the reduced matrix elements M^^^^^^^^ 



hard vertex function S^' 



IS 



{0*, V^* 



a 



(0) 



cos (4(y9, 



the double spin-flip 



(35) 



where 



and 
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1+1 
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entry 



(36) 
(37) 

proportional to 
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totic piece. 



b) produces an NNLO term in the unpolarized gg asymp- 



S''(^*,¥'*) 



a. 



Kl9 ® fa/hi 



(xi,/iF) 



Pg/g ® fg/h2 



{x2,fiF)- 



(38) 



The analogous quark function V^./i^^x, kx) corresponds to the transversity distribution 
and is odd under time reversal (T-odd). It cannot be generated radiatively through 
conventional PDF evolution from the T-even unpolarized function Vq^/h{x,kT) and does 
not contribute to the NLO asymptotic term. We flnd S' = 0, because the non-vanishing 



amplitudes Ti.q{qi^ 
(Ai = — A2). Therefore 



must have opposite helicities of the quark and antiquark 
the functions Vqjf^{x^h) contribute in pairs through the term pro- 
portional to S^'((y9*) = —{a'^e\'n/{2NcQ'^)) cos 293^. These contributions are anticipated to be 
much smaller than the usual spin-average contribution and negligible at large Q, in analogy 



to unpolarized Drell-Yan production |40l. l4l|. 



In summary, the azimuthal angle ((^*) dependence of photons in the gg scattering chan- 
nel is affected by large QCD contributions associated with interference between gluons of 
opposite helicities. These logarithmic corrections may arise at NLO through QCD radia- 
tion from conventional unpolarized PDF's, a mechanism that is unique to gluon scattering. 
Other types of spin-interference contributions (not considered here) involve spin-flip PDF's 
only. The soft and collinear logarithms associated with the spin-flip contributions must 
be resummed along the lines discussed in Ref. [Hi. Given that gg 77 is the sublead- 
ing production channel at the Tevatron^nd at the LHC, we henceforth neglect the gluon 
spin-flip contributions to the resummed W{Q^ b, y, fl^,), while subtracting the corresponding 
V5*-dependent asymptotic contribution from the flnite-order 2^3 cross section. The nature 
of gg spin-flip contributions can be explored by measuring the double-differential distribu- 
tion in (p^: and Qt at the LHC, a topic that is interesting also from the point of view of the 
Higgs boson search. Full resummation of the gluon spin-flip contributions may be needed in 
the future. 
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Figure 4: Comparison of [Pg/g'^ fg/p\{x, hf) and [P'g/g'^ fg/p]{x, hf) for the gluon PDF fg/p{x,HF) 
in the proton (multiplied by x^'^ to better illustrate the small-x region) at several values of the 
factorization scale ij,f- 



D. Complete expressions for resummed cross sections 

In this Section, we review complete expressions for the unpolarized resummed cross sec- 
tions, starting from the perturbative QCD approximation Wpert{Q,b,y,^*) valid at small 
impact parameters 6^ ^ 1 GeV~^. For a hard-scattering function J2hei.\'HiQj^*)f = 
^a{d*)h1{Q,9^), the form factor Wpert{Q,b,y,VL^) is 



X [Ca/ci O /ci/hi] (Xi, 6; /i) [C-a/c2 ® /ca/ha] {^2, b] fl) 



(39) 

The Sudakov function is defined in Eq. ( l24ll . and the function ha{Q,0^) collects radiative 
contributions to H{Q,9^) arising at NLO and beyond. We compute the functions ha, Aa, 

a ^, a nd Ug, respectively. The Aa, Ba and Ca/c coefficients 
45 . 66,[63l and listed in a consistent notation in Ref. Q. 

to join the small-Qr resummed cross sections W 
P — Ais the difference between 



Ba and Ca/c up to orders a„ 
are taken from Refs. 0, IH, \M 

We use a procedure outlined in Ref. [4 
with the large-Qr NLO cross sections P. In Eq. (|2T1) . Y 
the perturbative cross section P and its small-Qr asymptotic expansion A, explicitly given 
in Eq. f flGll . For each value of Q and y of the 77 pair, W + Y approaches P from above 
and eventually becomes smaller than P as Qt increases. We use W + Y as our prediction 
at Qt values below this point of crossing and the finite-order cross section P at Qt above 
the crossing point. 

The final cross sections depend on several factorization scales: Ci/6, C2Q, fi = C^/b in 
the W term, and fiF = C4Q in the Y term. Here Q {i = 1, ..4) are dimensionless constants 
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of order unity, chosen as C2 = C4 = 1, Ci = C3 = 2e~'^'^ = 1.123... by default. These choices 
simplify perturbative coefficients by eliminating scale-dependent logarithmic terms, cf. the 
appendix in Ref. Q. Dependence on the scale choice is studied in Section HVl 

In the general formulation of CSS resummation presented in 2^, 27 1, one has the freedom 
to choose different "resummation schemes", resulting effectively in variations in the form of 
ha{Q, 0^). These differences are compensated, up to higher-order corrections, by adjustments 
in the functions B and C. 

In "the CSS resummation scheme" one chooses ha{Q,0*) = 1, while including the 
virtual corrections to the 2-^2 scattering process in B and C. In this scheme, some B and 
C coefficients depend on the 2 — 2 hard scattering process and also on 6**. 

In an alternative prescription by Catani, de Florian and Grazzini ji^l, "the CFG resum- 
mation scheme", one keeps the 2 — 2 virtual corrections within a single function \H{Q, 9^:)f . 
In this case, the B and C functions depend only on the initial state. Most of our numerical 
calculations are realized in the CSS resummation scheme, with a few made in the CFG 
scheme for comparison purposes. 

In impact parameter (b) space used in the resummation procedure, we must integrate 
into the nonperturbative region of large 6, cf. Eq. (l22l) . Contributions from this region are 



known to be suppressed at high energies |48|, but some residual dependence may remain. 



In the qa + qg — > 77 channel, our model for the nonperturbative contributions (denoted 
as KNl [49|) is derived from the analysis of Drell-Yan pair and Z boson production. The 
nonperturbative function in this model is dominated at large Q by a soft contribution, which 
does not depend on the flavor of initial-state light quarks. This function is therefore expected 
to be applicable to the qq + qg ^ 77 process. 

The nonperturbative function in the gg + gqs channel, which is yet to be measured di- 
rectly, is approximated by the nonperturbative function for the qq + gq channel multiplied 
by the ratio Ca/Cf = 9/4 of the color factors Ca and Cp for the leading soft contributions 
in the gg and qq channels. This ansatz suggests stronger dependence of the gg + gqs channel 
on the nonperturbative input compared to the qq + qg channels. It leads to small differences 
from the prescription used in Refs. 0, 2^, where only the leading InQ term of the nonper- 
turbative function was rescaled. To examine the dependence of the resummed cross sections 
on the nonperturbative model, we evaluate some of them assuming an alternative (BLNY) 
parameterization of the nonperturbative function [s^l- 



IV. NUMERICAL RESULTS 



The analytical results of Sec. Ill are implemented in our computer codes Legacy and 
ResBos 46, 50, 51, 52|. We use the same parameters as in the calculation of Ref. and 
we concentrate on the region Qt < Q where our calculation is most reliable [5|. 



A. Results for Run 2 at the Tevatron 

In this section, we present our results for the Tevatron pp collider at = 1.96 TeV. 
We make the same restrictions on the flnal-state photons as those used in the experimental 
measurement by the Collider Detector at Fermilab (CDF) collaboration [Hll: transverse 
momentum > Prmin = 14 (13) GeV for the harder (softer) photon, and rapidity \y^\ < 0.9 
for each photon. We impose photon isolation by requiring the hadronic transverse energy 
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pp^yyX, VS = 1.96 TeV 
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Figure 5: Parton flavor decomposition of the resummed transverse momentum distribution at the 
energy \/5 = 1.96 TeV of the Tevatron Run-2. The total (solid), qq + qg (dashes), and gg + gqs 
(dash-dots) initial-state contributions are shown separately. 

not to exceed 1 GeV in the cone AR = 0.4 around each photon, as specified in the CDF 
publication. We also require the angular separation AR^^ between the photons to be larger 
than 0.3. 

We focus in this paper on the role of the gg contribution, referring to our other papers 
[sl, @| for a more complete treatment. 

To illustrate the relative importance of the individual initial-state contributions in the 
final answer, we provide a parton fiavor decomposition of our resummed transverse mo- 
mentum distribution da/dQT in Fig. El This distribution is integrated over all diphoton 
invariant masses Q, subject to the CDF cuts, and receives dominant contributions from the 
Qt < Q region. The gg + gqs contribution supplies about one-third of the total rate near 
Qt = 5 GeV. It falls steeply after Qt > 20 GeV, because the gluon PDF falls steeply with 
parton fractional momentum x. 

Dependence of the resummed cross sections on the choice of factorization scales men- 
tioned in Section IIII Dl is examined in Fig. [6l We pick a few characteristic combinations 
of alternative scales to probe the scale dependence associated with the resummed Sudakov 
function e""^, the 6-dependent PDF's Va/h{x,b) ^ [Ca/c ® fc/h]{x,b]fi), and the regular Y 
term. The small-Qr region is sensitive primarily to the scales Ci/6, C2 Q, C3/6 in the re- 
summed term W. The event rate at large Qt is controlled by the choice of the factorization 
scale = C4Q in the regular term Y. 

At the relatively low values of Q relevant for the Tevatron experiments, the scale depen- 
dence of the next-to-leading order gg + gqs cross section is still substantial, with variations 
being about -20% (+50%) at Qt = 5 - 10 GeV, ±10% at Qt = 10 - 20 GeV, and ±20% 
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Figure 6: Scale dependence of the cross sections in gg + gqs and all scattering channels. The upper 
frame shows cross sections for the default choice of scales specified in Section [ill Dl (solid), as well 
as for varied scales C2 = C4 = 2 (dashes), C2 = C4 = 0.5 (dots), and C3 = Ae~'^'^^ (dot-dashes). 
The lower two frames show ratios of the cross sections computed for the varied factorization scales 
to the cross section for the default choice of the scales. 



at Qt = 20 — 40 GeV. Since the Y term is the lovi^est-order approximation for gg 77c/ 
at Qt ~ Q, the scale dependence associated with the constant C4 remains pronounced at 
large Qt- The inclusive gg + gqs rate, integrated over Qt, varies by 20 — 40% almost in- 
dependently of the 77 invariant mass Q. The large scale dependence of the NLO gg + gqs 
cross section reflects slow perturbative convergence in gluon gluon scattering, observed also 
in other similar processes, e.g., gg —>■ Higgs via the top quark loop H, 55, H^. For this rea- 
son, a NNLO calculation would be desirable to reduce the scale uncertainty in the gg + gqs 
channel. 

On the other hand, the scale dependence of the cross section when all channels are 
combined is relatively mild, with variations not exceeding 10% at small Qt and 20% at 
large Qt- Variations in the integrated inclusive rate for all channels combined are below 
10% at Q > 30 GeV. 

Another aspect of scale dependence is associated with the assumed arrangement of log- 
arithmic terms in the resummed W term, i.e., the "resummation scheme" that is adopted. 
This dependence is yet another indicator of the size of higher-order corrections not included 
in the present analysis. Figure [7](a) shows ratios of the full resummed cross sections in the 
Catani-de Florian-Grazzini (CFG) and Collins-Soper-Sterman (CSS) resummation schemes. 
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Figure 7: Ratios of resummed cross sections at the Tevatron Run-2 computed in (a) the Catani-de 
Florian-Grazzini (CFG) and CoUins-Soper-Sterman (CSS) resummation schemes and (b) using the 
BLNY and KNl nonperturbative models, as functions of the 77 transverse momentum Qt- The 
ratios are shown in the qq + qg (dashed), gg + gqs (dot-dashed), and all (solid) scattering channels. 
A Qt < Q cut is imposed in this comparison. 



as described in Sec. III. The differences between these schemes stem from the different 
treatment of the NLO hard-vertex correction h^^\9^). The magnitude of /i^^)(6'^,) determines 
whether the channel is sensitive to the choice of the two resummation schemes. The mag- 
nitude of hSg\0^) in the gg + gqs channel exceeds that of h'^q\Oj^ in the qq + qg channel by 
roughly an order of magnitude for most values of the 9^ angle |43||. Consequently, while the 
dependence on the resummation scheme is practically negligible in the dominant qq + qg 
channel (dashed line), it can reach 15% in the subleading gg + gqs channel (dot-dashed 
line). The Qt spectrum in gg + gqs channel is slightly softer in the CFG scheme up to 
the point of switching to the fixed-order cross section at Qt ~ 60 GeV. The resummation 
scheme dependence in all channels (solid line) is less than 3-4%, reflecting mostly the scheme 
dependence in the gg + gqs channel. 
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To examine the sensitivity of the resummed predictions to long-distance nonperturbative 
dynamics in hadron-hadron scattering, we include in Fig. [7](b) a comparison with the re- 
summed cross sections for an alternative choice of the nonperturbative model. As explained 



in Sec. IIII Dl our default calculation is performed in the recent KNl model |49] for the 
nonperturbative part of the resummed form factor W{Q,b,y,Q^). Figure [7](b) shows ratios 
of the predictions for a different BLNY model [H^l and our default KNl model in various 
initial-state scattering channels. 

The difference is maximal at the lowest Qt, as expected, and it is less than 5% for 
the total cross section. For the qq + qg and gg + gqs initial states the maximal difference 
is about 5% and 20%, respectively. The dependence on the nonperturbative function is 
stronger in the gg + gqs channel, where the BLNY/KNl ratio in the gg + gqs channel 
reaches its maximum of 1.15 at Qt ~ 25 GeV and slowly decreases toward 1, reached at 
the switching point at Qt ~ 60 GeV. This behavior reflects our assumption of a larger 
magnitude of the nonperturbative function in the gg + gqs channel, which is rescaled in our 
model by Ca/Cf = 9/4 compared to the nonperturbative function in the qq + qg channel. In 
summary, despite a few-percent uncertainty associated with the nonperturbative function 
in the gg + gqs process, the overall dependence of the Tevatron 77 cross section on the 
nonperturbative input can be neglected. 

B. Results for the LHC 

To obtain predictions for pp collisions at the LHC at = 14 TeV, we employ the cuts 
on the individual photons used by the ATLAS collaboration in their simulations of Higgs 
boson decay, /i — >^ 77 We require transverse momentum pj^ > 40 (25) GeV for the harder 
(softer) photon, and rapidity \y'^\ < 2.5 for each photon. We impose the ATLAS isolation 
criteria, looser than for the Tevatron study, requiring less than 15 GeV of hadronic and 
extra electromagnetic transverse energy inside a AR = 0.4 cone around each photon. We 
also require the separation AR^^ between the two isolated photons to be above 0.4. The cuts 
optimized for the Higgs boson search may require adjustments in order to test perturbative 
QCD predictions in the full 77 invariant mass range accessible at the LHC. 

Distributions in the invariant mass Q, transverse momentum Qt, and azimuthal angle 
separation A(p = (p^-^ — tp^^ between the two photons in the laboratory frame are shown in 
Fig. m As before, we compare the magnitudes of the qq + qg and gg + gqs cross sections. 
The qualitative features are similar to those at the Tevatron, but the relative contribution 
of the various initial states changes at the LHC. The gg + gqs initial state contributes about 
25% of the total rate at Q ~ 80 GeV where the mass distribution peaks, but the gg + gqs 
rate falls faster than qq + qg with increasing invariant mass. 

In the invariant mass range relevant for the Higgs boson search, 115 < Q < 140 GeV, the 
transverse momentum distribution in Fig.[8](b) shows that the gg + gqs initial state accounts 
for about 25% of the rate at low Qt- At high transverse momentum, on the other hand, the 
other channels dominate. The relative size of the gg + gqs contribution drops as the invariant 
mass or the transverse momentum of the photon pair grows. The gg + gqs contribution falls 
more steeply with Qt for larger masses of the diphoton. These features are attributable to 
the steeply falling gluon distribution as a function of increasing momentum fraction x. 

The scale dependence at the LHC, presented in Fig. [9l is somewhat reduced compared to 
the Tevatron (cf . Fig. [6|) . Maximum scale variations of about 40% in the gg + gqs channel 
are observed at the peak of the da/dQT distribution, and they are substantially smaller at 
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Figure 8: Resummed da/dQ, da/dQr, and da/dA(p distributions of photon pairs at the LHC for 
ATLAS kinematic cuts. 



large Qt- The scale variation in the sum over all channels does not exceed 10% (15%) at 
small Qt (large Qt)- Variations in the integrated inclusive rate at Q > 50 GeV are below 
7% (30%) in all channels {gg + gqs channel). 

The dependence on the resummation scheme is mild at the LHC (cf. Fig. [TOlf a)). with 
the maximal differences between the CSS and CFG schemes below 0.5%, 10%, and 2% in 
qq+qg, gg+gqs, and all channels. The scheme dependence is again the largest in the gg+gqs 
channel, where it persists up to the point of switching to the fixed-order cross section at 
Qt ~ 120 GeV. The ratios of the resummed cross sections calculated in the BLNY and KNl 
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Figure 9: Scale dependence in the gg + gqs and all scattering channels at the LHC for the same 
scale choices as in Fig. [6l 



models for nonperturbative contributions in the CSS scheme are shown in Fig. [TOlfb). The 
influence of the long-distance (large-6) contributions is suppressed at the high center-of-mass 
energy of the LHC. Differences between the predictions in the two models do not exceed 
2%, 6%, and 2% in the qq + qg, gg + gqs, and all scattering channels. 

The KNl and BLNY nonperturbative models neglect the possibility of a strong x de- 
pendence of the nonperturbative function, which may substantially modify our predictions 
at the energy of the LHC collider. Analysis of small-x semi-inclusive deep inelastic scatter- 
ing data [57| suggests that x-dependent nonperturbative corrections of uncertain magnitude 
may substantially affect the resummed cross sections. Such corrections can be constrained 
by studying the rapidity and energy dependence of the nonperturbative function at the 
Tevatron and LHC, for example, from copious production of Z bosons [57]. We conclude 
that uncertainties due to the choice of the resummation scheme and the nonperturbative 
model will be small at the LHC, if the resummed nonperturbative function does not vary 
strongly with x. 



C. The role of the gqs contribution 

Figures ISlfTOl show the contributions from the qq + qg and gg + gqs channels along with 
their sum. One may wonder if a further decomposition into qq and qg (or gg and gqs) contri- 
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Figure 10: Same as Fig. (Tj at the LHC. 



butions could provide additional insights into the relative importance of different scattering 
processes. We observe in our calculations that the resummed cross sections W + Y and the 
fixed-order cross sections P in the elementary scattering subchannels (gg, qg ■,■■■) may not 
cross until Qt is significantly larger than Q. This result is at variance with our expectation 
that the fixed-order answer should be adequate when Qt is of order Q, where logarithmic 
effects are small, and the one-scale nature of the dynamics seems apparent. 

Consider, for example, the gg and gg -\- gqs transverse momentum distributions in the 
mass interval 115 < Q < 140 GeV at the LHC shown in Figs. [TTlfa) and (b). In the gg 
channel alone (Fig. [TTVa)). the W + Y cross section remains above the NLO cross section P 
until Qt ~ 140 GeV. However, after the gqs contribution is included (Fig. [TlT b)). W + Y 
crosses P at Qt ~ 105 GeV. Our expectation of the adequacy of the NLO prediction at 
Qt ~ Q is satisfied in this case, and this conclusion also holds for other intervals ofQ. At the 
crossing point, the two cross sections satisfy W + Y = P, i.e., W = A; the resummed term 
is equal to its NLO perturbative expansion, the asymptotic term. Similarly, good matching 
of the resummed and NLO cross sections in the qq + qg channel requires that we include 
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Figure 11: Inclusion of the qg contribution improves the matching of the resummed and NLO 
perturbative cross sections at large Qt, as demonstrated by these plots of the resummed and finite- 
order NLO cross sections for (a) the gg channel only; (b) the combined gg + gqs channel. The 
resummed and NLO cross sections are shown by the solid and dashed lines. 

both qq and qg contributions. 

This feature can be understood by noticing that the flavors of the PDF's fa/h{x,fi) 
mix in the process of PDF evolution. Consequently the perturbative expansion of W 
in the gg channel contains the full NLO asymptotic piece A in the combined gg + gqs 
channel, generated from the Sudakov exponential and lowest-order resummed contribution 
oc fg/hi{xi, 1/6) fg/h2{x2, 1/6) evaluated at a scale of order 1/6. The mismatch between the 
flavor content in the perturbatively expanded W and A in nominally the same gg subchannel 
causes the difference 1^ — A to be large and delays the crossing. On the other hand, the 
flavor content of W and A is the same (up to NNLO) when the gg and gqs contributions 
are combined, and the matching is improved. The gqs scattering subchannel has been as- 
sumed to be small and neglected in past studies, and indeed it contributes about one tenth 
of the gg + gqs inclusive rate da/dQ. However, we see that the gqs contribution must be 
included to correctly predict da/dQx and to realize matching between the resummed and 
perturbative contributions at large transverse momenta. 

V. SUMMARY AND CONCLUSIONS 

In this paper, we address new theoretical issues in Qt resummation at two-loop accuracy 
that arise in the gluon-gluon subprocess, gg + gq ^ 'J'J, one of the important short-distance 
subprocesses that contribute to the inclusive reactions pp —>■ 77X at the Fermilab Tevatron 
and pp 77X at the CERN Large Hadron Collider (LHC). 

We evaluate all next-to-leading (NLO) contributions of order 0{a'^al) to the gg+gq — > 77 
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process (Fig. [I](b-e)). A new ingredient in this paper is the inclusion of the gq 775 
process, Fig. [T](d), a necessary component of the resummed NLO contribution. We resum 
to next-to-next-to-leading logarithmic (NNLL) accuracy the large logarithmic terms of the 
form ln(Q|./(5^) in the limit when Qt of the 77 pair is smaller than its invariant mass Q. 
The perturbative Sudakov functions A and B and the Wilson coefficient functions C in the 
resummed cross section W are computed to orders a^, a^, and as- The resummed cross 



sections are computed according to the CSS [1| and CFG [43| resummation schemes, with 
the differences between the two approaches reflecting the size of higher-order corrections. A 
new nonperturbative function [4^, dominated by a process-independent soft correction, is 
employed to describe the dynamics at large impact parameters. 

Subtraction of the singular logarithmic contributions associated with initial-state radia- 
tion from the NLO cross section P defines a regular piece Y. This regular term is added to 
the small-Qr resummed cross section W to predict the production rate at small to moderate 
values of Qt- In the gg channel, we also subtract from P a new singular spin-fiip contribu- 
tion that affects azimuthal angle dependence in the Collins-Soper reference frame. For 
our final prediction, we switch from the resummed cross section W + Y to P at the point 
where W + Y crosses P, approaching P from above, as in Ref. [i^. The location of this 
point in Qt is of order Q in the qq + qg and the gg + gq channels. For such matching to 
happen, it is essential to combine cross sections in the qq and qg {gg and gq) channels, as 
demonstrated in Sec. IIV Ci 

At the LHC (Tevatron), the gg + gq subprocess contributes 20% (10%) of the total 77 
production rate (integrated over the full range of the photons' momenta). The relative 
contribution of gg + gq scattering may reach 25% for some Q and Qt values. The gg + gq 
channel provides an interesting opportunity to test CSS resummation at a loop level and may 
be explored in detail at later stages of the LHC operation. The NNLL/NLO resummed cross 
section for the gg + gqs channel is used in Ref. ^ to predict fully differential distributions 
of Higgs bosons and QCD background at the LHC in the Higgs — > 77 decay mode. 
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Appendix A: THE gqs 'j-fqs AMPLITUDE 



To obtain the gluon-quark contribution to the gg + gqs scattering channel shown in 
Fig. [11(d), we derive the helicity-dependent qqjjg amplitude A^5(i?i, ^2, 73, 74, S's) from 
the one-loop qqggg amplitude in the color-decomposed representation available in Ref. 
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The qq'yyg amplitude is expressed as 



MM\Q2\l3,l^,9r) = ^9'e'{J24\T^,% E 4'^''(l,,2,-;a(3),a(4),a(5)) (Al) 



in terms of the primitive amplitudes A^.'l^^'^\lq,2q]3,4,5) = —Al.i{lq,2q; 3,4,5) — 
^5.]^(lg, 2g; 3, 4, 5) for qqggg — > scattering involving a spin-1/2 fermion loop. The am- 
plitude A^5 is proportional to the sum X^i; (^^^^J of squared quark charges circulating in 
the fermion loop, as well as the QCD generator matrix T^^^^, with Tr(T"iT"2) = 5°'^°''^ . The 
color indices Ci, C2, and belong to the quark 1, antiquark 2, and gluon 5. The primitive 
amplitudes are summed over all possible permutations Sf'^^^ of the legs 3, 4, and 5. 

Equation ( 1A1[ ) is derived from Eq. (2.10) of Ref. Q after gluons 3 and 4 are replaced 
with photons, i.e., the QCD generators T°^^djvd T"* are replaced by identity matrices, and 
the overall charge factor is adjusted, g^ '^g^J^iii^^n)^- It correctly reproduces the small- 
Qt asymptotic behavior reflected in Eq. (ITGll . which we derive by applying factorization 
relations in the splitting amplitude formalism discussed in Appendix [Bl The qq'yyg am- 
plitude in Eq. (lAip disagrees with the one published in Ref. 2J| which appears to violate 



factorization relations in the q \\ q limit. A few independent ainplitudes ^{.^(Ig, 2^; 3, 4, 5) 
and ^^.^(Ig, 2^; 3, 4, 5) are presented explicitly in Sec. 5 of Ref. |8|, with the remaining am- 
plitudes related by discrete symmetries according to Eq. (5.25) in that publication. Some 
qqggg amplitudes contain infrared poles, which cancel in the sum over permutations S^^'^^\ 
We retain only non- vanishing finite parts of such divergent amplitudes, i.e., we take 
Af.^ = iF^/ (167r^) for x = / and s. 



Appendix B: DERIVATION OF THE SMALL-Qt ASYMPTOTIC TERM FOR 
GLUON-GLUON SCATTERING 

In this appendix, we derive the small-Qr asymptotic approximation Eq. (fT6l l for the 
NLO cross section in gig2 7374 scattering. We expand the finite-order cross section as a 
series in the small parameter Q'^/Q'^. Consider first the leading real-emission contributions, 
which arise when gluon 5 is radiated off the external gluon leg 1 or 2 as in Fig. [I](b).^ In 
the notation introduced in Sec. Ill, the small-Q^ approximation for the real-emission cross 
section takes the form 

/■I /-i 1 1 2 

A{Q,QT,y,n*)l^^l = j^^dii J^^d^2f9/h^{^l,^^F)fg/h2{^2,^iF)J^-^^^;^\M5\ 

The right-hand side of Eq. (IBlll includes a product of the gluon parton densities 
fg/hi{^i,2, fJ'p), squares of parton-scattering amplitudes A^s, and phase-space factors, in- 
tegrated over the light-cone momentum fractions ^1,2 = P^2/-fV2 of the incoming gluons 1 



^ In contrast, Feynman graphs with gluon radiation off a propagator in the quark loop [Fig.[ljc)] are finite 
in the Qt limit. 
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and 2. The delta- functions constrain integration to phase-space regions where the final-state 
gluon 5 is collinear to gluon 1 [^5 = (1 — Xi)pi], coUinear to gluon 2 [ps = (1 — 22)^2)5 or 
soft 0], with Xi = Xi/^i for i = 1,2. 

The 2 — > 3 helicity amplitude A^5(l, 2, 3, 4, 5) is analyzed conveniently in an unphysical 
scattering channel g{Pi,^i)giP2,^2)l{P3,\)l{P4,^4)g{P5,^5)- The momenta and 
helicities Aj are related to the physical momenta and helicities Aj as {pi, Xi}={—pi, —Xi} 
for i = 1 or 2, and {pi, Xi}={pi, Xi} for i =3, 4, or 5. A^5(l, 2, 3, 4, 5) is a shorthand notation 
for A45(pi, Ai;p2, A2;P3, A3;p4, A4;p5, A5). 



The amplitude A^5(l, 2, 3, 4, 5) was derived in Refs. 0, mi from color-decomposed 5- 
gluon 1-loop scattering amplitudes Al5(l, 2, 3, 4, 5) is built from 1-loop partial ampli- 
tudes y45;i(l, 2, 3, 4, 5) for the ggiig scattering process, identical to the partial ampli- 
tudes for ggggg scattering via a spin-1/2 fermion loop [9|. The squared 5-leg amplitude, 
averaged over spins, colors, and identical final-state particles, is 



a 



(1) 



_ _E_ _ 

'^li-'l2,A3,A4,A5 



^ 0'2, (73, Cr4, (T5) 



(B2) 



with 



(1) 



2\2 



Nr. 



m - 1' 



(B3) 



The partial amplitudes are summed over all permutations a of the external indices (1,2,3,5) 
with a fixed cyclic ordering of (1,2,5), i.e., cyclically-ordered (COP) permutations: 

^ A;l(0"l, 0"2, Cr3, (T4, 0-5) = 

A5;i(l, 2, 5, 3, 4) + ^5.1(1, 2, 3, 5, 4) + A,,,{1, 3, 2, 5, 4) + A,,,{3, 1, 2, 5, 4) + 
A5;i(5, 1, 2, 3, 4) + A5;i(5, 1, 3, 2, 4) + A,,,{5, 3, 1, 2, 4) + A^.i^S, 5, 1, 2, 4) + 
A5;i(2, 5, 1, 3, 4) + A5;i(2, 5, 3, 1, 4) + A,,^i2, 3, 5, 1, 4) + A,,,{3, 2, 5, 1, 4). (B4) 

The collinear and soft behaviors of the amplitude can be established by followin g th e 
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approach in Refs. extended recently to the two-loop level [5 

When gluon 5 is collinear to gluon 1, the amplitude A^5(l, 2, 3, 4, 5) is dominated by six par- 
tial amplitudes with cyclically adjacent indices 5 and 1, such as A^-i^b, 1,2,3,4). Similarly, 
when gluon 5 is collinear to gluon 2, Ad^^l, 2, 3, 4, 5) is dominated by six partial amplitudes 
with cyclically adjacent indices 2 and 5. Each leading partial amplitude A^.i{...,5,l, ...) 
factors in the 5 || 1 collinear limit into a 4-leg partial amplitude ^44.1 /, . . . ) for production 



of 2, 3, 4, and intermediate gluon J, and amplitude Split 
tree-level splitting of / into 5 and 1: 



tree 
A 



!(5 



1) (mI, S [e; 



describing 



A 



5,1, 



5 1 



A7=±l 



spiit*:f r5, 1) A 



4:11 



/, ...) + subleading terms. 



(B5) 



The ellipses in Eq. f lB5l ) denote the same permutation of indices 2, 3, and 4 in Ar,.i and A^-i. 
The amplitudes Split^|^(5, 1) are universal functions of the momenta pi, pi, and ps, which 
in our case satisfy pj = pi + p^, pi = {1 — z)pi, and p^ = zpj, where z = 1 — l/xi. The 
right-hand side of Eq. (|B5l) is summed over the helicities Xj of /. The collinear factorization 
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relation applies to any one-loop n— leg primitive amplitude A^°°p{1, ... 



n): 



SplitiT,(«,&)^n-"i(---,/,--0 + SplitLl^ (a, /,...) 



(B6) 



Equation f lB6p is evaluated here for = 5 external legs, along with the condition that the 
tree primitive amplitude Af'^^ vanishes in — >^ 9911 process. 

Using Eqs. (1B2I) . (IB4I) . (1B5I) . we derive the approximate form for I^Wsp in the 5 || 1 limit: 

|Al5(l, 2, 3, 4, 5) I' ^ a« J] A^Kp,, A^; 2, 3, 4)T^, ,;,^(xi)^4(p/, A,; 2, 3, 4). (B7) 

Here A^4(/, 2, 3, 4) = Z^crg^g A^-i^I, a2, crs, <J4,) is the normalized 4-leg amplitude, obtained by 
summation of the partial amplitudes ^4.1(1, 2, 3, 4) over all possible permutations 5*3 of the 
legs 2, 3, and 4. The amplitude A^4 and complex-conjugate amplitude A41 are evaluated 
for independent helicities A/ and \'j of I. Ty^ Xii^i) absorbs contributions from the splitting 
amplitudes: 



Ai,A5=±l 



Spiit*:f.(^;i,5) 



Spiit*_i;,(^;i,5). 



(B8) 



In a basis with A = +1 and A = —1, Ty^ Xj{x) is a matrix of the form 



^A.,A^(^) 



2C, 



2xpi ■ P5 



-^ + ^ + xn-x) -— M 

l—x X ^ ' X (51) 

-— w j^ + ^ + xn-x) 

X [51\ l—x X ^ ' 



(B9) 



The diagonal entries of T\^ yi{^\) give rise to terms proportional to the unpolarized splitting 
function Pgjgi^x) in the asymptotic cross section, with 



P,/3(xi) = 2Ca 



X\ 



\ —X\ 



+ - Xx, 



+ 



6 



<5(l-xi), (BIO) 



where Nf is the number of active quark flavors. The off-diagonal entries give rise to terms 
proportional to the spin-flip splitting function 



P^/^(xi) = 2C^(l-xi)/xi, 



(Bll) 



multiplied by the ratio of spinor products (51) = (5 + |1— ) and [51] = (5 — |1+). In a 
general reference frame, (51)/[51] is a complex phase depending on the azimuthal separation 
V^i — V's between the gluons 1 and 5. In the Collins-Soper frame, this phase reduces to 
(51)/[51] = -1.2 



A collinear approximation for jA^sl is derived in Ref. [23| in the framework of the dipole factorization 



formahsm 



25|. This approximation agrees with ours up to phases of the off-diagonal terms, which are not 



the same as in Eq. (|B9p . Our expression is shown upon a closer examination to produce correct phases in 
an arbitrary reference frame jg^. 
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Next, we employ explicit expressions for A^4(/,2,3,4) from Ref. [22|, given by products 
A44(/,2,3,4) = 'S'x^ A2 A3 A4 ^1^X2X3X4 reduced matrix elements m|^^^_^^^^^ and phase factors 
SxjX2X3>^4- With these expressions inserted, Eq. ( 1B7I1 becomes in the Collins-Soper frame 



1-^5(1,2,3,4,5)1 
where 



2xipi ■ P5 



LM) = _ E _ 

Ai,A2,A3,A4 



M- 



(1) 



A7A2A3A4 



and 



L'^ie,) = -4Re {Mi]l_,^^, + M«i,i,_i + M%^,^_, + l} 



(B12) 

(B13) 
(B14) 



In the 5 II 2 collinear limit |A^5|^ is 



|A^5(1,2,3,4,5)|2-^^^-^— |p,/,(x2)L,(^,) + P;/^(x2)L;(^^,) cos2(^J. (B15) 

In the soft limit P5 0, lA^sl^ factors as 

\M5{1, 2, 3, 4, 5) 1^ ^ -^2C^— |;V/l4(l, 2, 3, 4) f + subleading terms. (B16) 

Qt ^ 

Inserting collinear and soft approximations (1B12I) , (IB15I1 , and (1B16I) in Eq. (IBll l and making 
simplifications, we derive the asymptotic expression for real-emission contributions, 



some 



da. 



gg 



dQ'^dy dQ'^dVt^ 



real 



1 as 
S 27rQ| TT 



X 



+ 



Pg/g ® /fl/^lj (3^1, /^f) fg/h2{X2, fJ'p) + fg/hi{Xl, fip) [Pg/g ® /<;/h2j {^2, fip) 

+ fg/hA^lif^p) [Pg/g^fg/hi] {x2,^^f) 



(B17) 



Once we add the two-loop 4-leg virtual corrections [Fig. [T](e)], the soft singularities in the 
real-emission cross section residing at Qt = are canceled [231, l43||. The final small-Qr 
expression coincides with Eq. ( flGll . 
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